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Abstract — We consider a Gaussian MISO wiretap channel, where a 
multi-antenna source communicates with a single-antenna destination in 
the presence of a single-antenna eavesdropper. The communication is 
assisted by multi-antenna helpers that act as jammers to the eavesdrop- 
per. Each helper independently transmits noise which lies in the null 
' space of the channel to the destination, thus creates no interference to 
the destination. Under the assumption that there is eavesdropper channel 
uncertainty, we derive the optimal covariance matrix for the source signal 
so that the secrecy rate is maximized subject to probabiUty of outage 
and power constraints. Assuming that the eavesdropper channels follow 

■ zero-mean Gaussian model with known covariances, we derive the outage 
probability in a closed form. Simulation results in support of the analysis 
are provided. 

' Index Terms — Gaussian MISO wiretap channel, cooperative jamming, 
artificial noise, outage probability 

I. Introduction 

■ Physical layer secrecy exploits channel conditions to maximize the 
rate of reliable information delivered to the legitimate destination, 

■ with the eavesdropper being kept as ignorant of that information 
as possible. This line of work was pioneered by Wyner |T|, who 

, showed that when an eavesdropper's channel is a degraded version 
of the source-destination channel, the source and destination can 
exchange secure messages perfectly at a non-zero rate, while the 
eavesdropper can learn almost nothing about the messages from its 
observations. The need for physical layer security in the context of 
wireless communications is motivated by challenges associated with 
classical cryptographic approaches, most notably the exchange and 

■ maintenance of private keys. 

The secrecy capacity for multiple antenna wiretap channels with 

■ perfect channel state information on the eavesdropper is established 
\ in ||2|, fSj, m under sum power constraints, and in {5j and ||6) 

■ under power-covariance constraints. For MISO wiretap channels, the 
optimal input covariance matrix that achieves the secrecy capacity 

, under a sum power constraint was given in a closed form in (7|. 
As shown in [H, when the channel to the destination encounters 
more fading than the channel to the eavesdropper, a positive secrecy 
rate can not be guaranteed. One way to overcome this problem is 
to use helpers who amplify-and-forward, or decode-and-forward the 
source signal, or perform cooperative jamming (CJ). In the latter 
case, helpers do not need to receive nor relay the source message, but 
rather just transmit noise to degrade the channel to the eavesdropper, 
therefore increasing the secrecy rate. In (9), a multiple access wiretap 
channel was considered and it was shown that if the optimal power 
allocation policy does not allow a certain user to transmit, that 
particular user could increase the secrecy rate by transmitting artificial 
noise. In 1101 . multiple helpers were employed to transmit a weighted 
jamming signal that enforces nulling at the legitimate destination 
and maximize the secrecy rate. Following that work, II II found the 
optimal weights by avoiding the nulling at the destination, which 
achieves higher secrecy rate. 

In 1121 . A MISO system with multi-antenna transmitter and single- 
antenna legitimate receiver and eavesdropper was studied. The trans- 



mitter constructs a Gaussian distributed artificial noise that lies in 
the null space of its channel to the legitimate receiver, and sends a 
sum of information bearing signal and the artificial noise, which only 
confuses the eavesdropper but not the legitimate receiver Since the 
transmitter does not know the channel to the eavesdropper, the power 
of the artificial noise is designed to uniformly spread along the null 
space. In |13| , the use of artificial interference for a MIMO wiretap 
channel is studied. The transmitter sends a sum of signal and artificial 
noise. When the CSI of the eavesdropper is unknown, the artificial 
noise is designed to lie in the null space of the right singular vector 
associated with the largest singular value of its transmitter-receiver 
channel matrix, and the power of the noise is uniformly spread along 
the null space. In 1141 . the artificial noise is studied in a MISO 
channel overheard by multiple single-antenna eavesdroppers. The 
source transmits a mixture of message and artificial noise. Without 
the CSI of the eavesdroppers, the outage based artificial noise design 
is formulated, with a so-called safe convex approximation is used to 
find the solution. 

In this paper we consider a MISO wiretap channel with multiple 
multi-antenna helpers implementing cooperative jamming. A multi- 
antenna source transmits the message, while each helper transmits 
jamming noise that lies in the subspace that is orthogonal to the 
helper-destination channel. Each helper generates jamming noise 
locally based on only its own link to the receiver Without any other 
information, the power of the jamming noise is uniformly spread 
along the null space. We study the problem of determining the input 
covariance matrix so that the secrecy rate is maximized subject to a 
sum power constraint and also an outage probability constraint. The 
source can determine the optimal input covariance matrix solely based 
on statistical information about the source-eavesdropper channel, with 
the channel following a zero mean Gaussian model with identity 
covariance. Assuming that the helper-eavesdropper channels follow 
zero-mean Gaussian distributions with arbitrary covariance matrices, 
the outage probability isit is as obtained in a closed form. Introducing 
the outage constraint not only provides quality control but also 
allows taking the uncertainty of the eavesdropper's channel into 
consideration and simplifies the optimization problem with respect 
to the covariance matrix of the input. 

Notation - Throughout this paper, following notations are adopted. 
Upper case and lower case bold symbols denote matrices and vectors, 
respectively. Superscripts *, T and f denote respectively conjugate, 
transposition and conjugate transposition. Tr(A) denotes the trace of 
the matrix A. A denotes that the matrix A is Hermitian positive 
semi-definite. \a\ denotes absolute value of the complex number a. 
||a|| = V at a denotes Euclidean norm of the vector a. I„ denotes 
the identity matrix of order n (the subscript is dropped when the 
dimension is obvious). C" denotes the set of all n x 1 complex 
vectors. E{-} denotes the expectation operator i = \/— T. a- ~ y 
denotes x and y have identical distributions. 
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Fig. 1. System model. 



II. System and Signal Models 

We consider a Gaussian MISO wiretap channel with A'^ helpers, 
shown in Fig. [T] The transmitter Alice uses Nt antennas to send 
messages to the legitimate receiver. Bob, through the channel G 
C^' . The eavesdropper. Eve, intercepts messages through the channel 



go G 



The transmitter is aided by A'^ helpers; each helper has Nk 



antennas {Nk > 2,fc = 1, 



, N) and transmits noise to confound 



Eve. The channel from helper k to Bob is denoted by G C , fc = 
1, ■ • • , TV, and the channel from helper k to Eve is denoted by g^ G 

c^^fc = I,-- - ,A-. 

At Bob or Eve, the received signal is a combination of the source 
signal, jamming noise and additive white Gaussian noise (AWGN). 
The received signal at Bob and Eve can be expressed as: 



Vb 



5Z '^I'^fe + '^i" 



^goX + l]gInfc 



(1) 



(2) 



where n^, k = I,-- - ,A'' represents the jamming noise that is 
generated by the helpers; and x is the A'^t x 1 source signal vector with 
input covariance matrix 0; Ub and rie are the AWGN received 



at Bob and Eve with Eflrib] 1 = 



] = A^o- The source power 



is constrained as Tr(Ra;) < Ps. We may represent R^; = P^Q with 
Q and Tr(Q) < 1. The power available at helper k is P^. 

In the following, we assume that each helper knows only its own 
link to the destination, hfe.fc = 1, ■ • ■ , A'^, and locally designs its 
jamming noise so that it delivers a null at Bob. 

A. Nulling Noise Structure 

In order for the noise vector of helper k to cause nulling at Bob, 
it should hold that 



0, k 



(3) 



It is clear that each helper should be equipped with A'^^ > 2 antennas 
such that there are enough degrees of freedom to design rife based 
on l[3j- The general solution of l[3j can be expressed as rifc — E^Vfe, 
where Efc (Nk x {Nk - 1)) is the null space of h[,, with E[,Efc = I 
and Vfc is any arbitrary (A'fc — 1) x 1 vector. Let the jamming noise 
be rifc — WfcEfctfc, where Wk is a weight that will be selected to meet 
the power constraint, and ~ CA/'(0, 1). The power of the jamming 
noise is 

E{||nfep} = E{TrKnt)} = \wkfE{T^i^ktktiEl)} (4) 
= !wfe|'Tr(EfeEl) = {Nk - l)\wkf . 

The power constraint is hence (A'^fc — 1)1™^^ < Pk. 



With this nulling noise, the received signal at Bob and Eve becomes 
respectively. 



ye 



/AhJjX + rib, 

JV 

/Ag?,x + Vg[,(™feEfetfc) + ; 



fe=i 



(5) 
(6) 



The secrecy rate of this system equals to 



Ci = log2 ( 1 + ^^hjQho 



log2 1 + 



ELikfcPI|ElgfcP + iVo 



(7) 



To make Ci as large as possible, the maximum relay power should 
be used, i.e., it should hold that {Nk - 1)1^^1^ = Pk. With this, we 
have 



Ci = log2 (^1 +puh5Qho 
- log2 f 1 + 



PogoQgo 



(8) 



ELiPfe(A^fc-l)-M|ElgfeP + l 

where po = Ps/No, pk ~ Pk/No, fc = 1, ■ • ■ , TV are the Signal-to- 
Noise Ratio (SNR) at Alice and the fcth helper. 

In the following section we study the problem of maximizing the 
secrecy rate with respect to Q using only statistical information about 
the eavesdropper channels. 

III. Outage Constrained Secrecy Rate Maximization 

Let us assume that the nodes do not have any channel information 
with the following exceptions: 

• Alice knows ho perfectly 

• Alice has statistical information on go, i.e., go ~ CAf{0, a^I). 

• Helper k knows its own link to Bob, hfc. 
The outage probability is defined as 1151 



Pout(-R) = min Pr(Ci < R). 

q^O, Tr(Q)<l 



(9) 



We will determine the maximum R such that the outage probability 
is below a prescribed small level, e, which represents the qualify of 
service (QoS). Mathematically, the problem is formulated as 



max R 
Q 

s.t. Pout (i?) < e. 



(10) 



The problem of l IlOl l is equivalent to the problem of maximizing the 
secrecy rate subject to QoS and power constraints as follows 



max R 

q 

s.t. Q ^ 0, Tr(Q) < 1, 
Pr(Ci <R)<e. 

A. Optimal Input Covariance Matrix Structure 



(11) 
(12) 



Lemma 1: For the problem of i T/Tt . and for go ~ CM{Q, cr^I), 
the optimal Q is given by Q* = hohg/||hop, and the optimization 
problem can be written as 



max R 

R 

S.t. Pr 



Polffoi 



ElLiPfe(^fc-i)-M|Elgfc 



(13) 

.>^+^"lfrll'-ll<.. 



2« 



(14) 



The proof is given in Appendix IaI 

Next, we solve tiie problem of il3l . Obviously, for the optimal R, 
the constraint of l ll4b holds with equality. Let us define the critical 
value Xe so that 

12 



Pr 



(15) 



^ELiP4A^fc-l)-M|ElgfcP + l 
Then, the optimal R is given by 

R* = log(l + pollholl') - log(l + xe). (16) 
The condition for i?* > is given in the following lemma. 
Lemma 2: For a given e, _R* > if and only if 



Pr 



po\goi\ 



>po||ho|r <e. (17) 



To find the critical value Xe, defined in l ll5t , we can use the bisection 
method. The calculation of the outage probability is discussed the 
following subsection. 

B. Closed Form Outage Probability 

Let us assume that go ~ CA/'(0,(j^I) and g^ ~ CA/'(0,Sfe). In 
this case, the probability of outage can be found in a closed form as 
follows. 

The calculation of the outage probability involves calculating 



Pr 



Po\goi\ 



ELiPfe(A'fe-l)-M|Elgfer + 1 



> X 



(18) 



which can be rewritten as 

P^fEL P^^^^ - l)"'l|Elgfef - -\90i\' < -l). (19) 
\ fe— i X ^ 

We use the methods in 1161 to calculate l |19b . Note that 
Etgfe ~ CAA(0,Ej:SfcEfc). Let E^S^Efc have eigen-decomposition 
JfeDfeUt. 



D = diag( 



Pi 



rDi, 



Pn 



' 'iVjv-1 
After a few derivations, l |19l l is equivalent to 

Pr(z+Dz < -1) 



-a'po/x). (20) 



(21) 



where z ~ CAf{0,lN+i ). Let Y = z^Dz. Y is an indefi- 
nite quadratic form. The results in 1161 give the expression for 
Pt{Y > y). Denote ~ f^Po/x- L^t vi, - ■ ■ ,vk denote different 
positive diagonal entries of D, with multiplicity mi, • • • ,mK and 
mi + • • • + rriK ~ N. According to the result in 1161 Eq. (32)] for 
the case z ~ CA/'(0, Ijv+i), we have 



Pr{Y>y) = -J2 



CXp(-y/i/fc) -(mfc-l) 



where 

gk{s,y) = e'"^ 



(-iyfc)'"fc-i(mfc - 1)! 



n 



5fc 



(1 + £.o/Vk + Cos) fj^ (Qfej - !/jS)'3'=3 



\o,y) (22) 



(23) 



Ofej = 



-1 j = k 
rrij j 
1 j^k ' 



(24) 
(25) 



and 5^'"''"^'' {s,y) denotes the [ink — 1) order derivative of gk{s,y) 
with respect to s. Using the recursive computation method in 1161 
Eq. (21)], we have 



-(") 
9 k 



with 



(...)-e("7' 

1=0 \ 



gl!\s,y)[lngk{s,y)]^"-'\ n> 1 (26) 



[^ngk{s,y)Y ' = -y - 



Co 



1 + Co/i^fe + Cos ' Qfej — UjS 



(27) 



(n-l)!^r 



for n > 2. (28) 

IV. Simulation AND Analysis 
In simulation, Alice has three antennas, each helper has two 
antennas. Alice-Eve link go ~ CA/'(0,I), and kth helper-Eve link 
gfe ~ CA/'(0,I). Alice has SNR po = 5 dB, and all iV helpers 
have SNR pk = 2 dB. Outage probability constraint, e = 0.01, and 
Pout{R*) ~ e, where R* is given in ( 116b . Results are averaged 
over 10^ independent trials. For each trial, generate independent and 
identical distributed ho, ho ~ CA/'(0, 1). During one trial, ho is kept 
as constants. For a R* > 0, lemma [2] has to be satisfied. When 
ho is much more degraded than go, with small number of helpers, 
the secrecy rate can not be guaranteed to be positive. During each 
trial, for = 5, • ■ • , 10 number of helpers, use bi-section method 
to search for a larger R that satisfying the outage constraint, until 
it converges to R* . In Fig. (2] the secrecy rate increases when the 
number of helpers increases. In Fig. [3] since Pr(i?) is a function of 
R, fix Ri = 0.61og2(l + Ps\\hof), plot the Pr(_Ri), the outage 
probability decreases when the number of helpers increases. Because 
Bob is not affected by the noise transmitted by the helpers, but Eve 
is, and the more helpers, the more confounded Eve will be. 
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Fig. 2. Secrecy Rate vs Number of Helpers, po = 5 dB, pk = 2 dB. 

V. Conclusions 

We proposed a cooperative jamming scheme, where multiple 
helpers transmit nulling noise to maximize secrecy rate subject to an 
outage probability constraint, and a power constraint. Assuming that 
the transmitter only knows its channel to legitimate receiver and the 
statistical CSI of its channel to eavesdropper, each helper known its 
own link to receiver, we formulated and solved an outage constrained 
secrecy rate maximization problem. Simulations show that proposed 
design could guarantee a low outage probability. 




Number of helpers 

Fig. 3 . Outage Probability vs Number of Helpers when fix i? = . 6 log2 ( 1 + 
Ps\\hof). po =5 dB, pk =2 dB. 



Appendix A 
Proof of Lemma[T] 



After a few derivations we get 



Pr(Ci < i?) = Pr 



POgpQgO 



ELiP^(^fe-l)-M|ElgfcP + l 



1 + pohtQho 

2« 



(29) 



To determine the structure of the optimal Q, let Q have an eigen- 
decomposition Q = UAU^ where A — diag(Ai, ■ ■ ■ , Ajvj) is the 
matrix of eigenvalues with Ai > ■ ■ ■ > Ajvj and U is an unitary 
matrix. First, we notice that 



g^Qgo = (Utgo)tA(Utgo) ~ gjAgo = ^"^^ A,|5o.|' (30) 

where we have used the fact that if go ~ CJ\f{0, o-qI), then U^go ~ 
go. According to l l29t and l l30t . we have 

pogjAgo 



Pr(Ci <R) = Pr(^ 



j:tip>=iNk-i)-'\\Eisk\\~ + i 

^ 1 + pohJUAUtho 



2« 



- 1 



(31) 



From d31t . we know that for the optimal A, the optimal U should 



maximize hJUAUt no, since 



1 is increasing 



with hQUAU^ho but decreasing with R, in other words, a larger 
hJUAU^ho will allow a larger R without violating the outage con- 
straint Pr(Ci < R) < e. From this fact, we can reveal the structure 
of U. Let Utho = y = [j/i,-- - ,yNtf. Then ||y||2 = ||hof . With 
this, we write 



hJUAUtho = ytAy < Ai||yf = Ai||ho| 



(32) 



Equality holds if y = [||ho||, 0, ■ • ■ , 0]^. It follows from U^ho = y 
that U = [ho/||ho||, U2, • • • , UNt]- According to lHOll, (ED and ( [32] l. 
we have 



Pr(Ci<i?) = Pr 



ELiPfe(A^fe-l)-M|ElgfeP + l 
1 + poAi||hof 



> 



2« 



(33) 



From ( 133b . we know that for the optimal Ai, smaller values of 
A2 , • • • , AjV( will allow a larger R without violating the outage 



constraint Pr(Ci < R) < e. Thus, it should hold that A2 = • ■ • = 
AjV( = 0. As a result, we have 

|2 

Pr(Ci<i?) = Pr( 



po\goi\ 



ELiP^(^^-l)"M|Elgfe|P + l 
Pollhof 1 



From i34l . we know that 



> 2« 

Poll hp f 



1 

2" 



Ai 



(34) 



is increasing with Ai but 



2« Ai 

decreasing with R, in other words, a larger Ai will allow a larger 
R without violating the outage constraint Pr(Ci < R) < e. Since 
Tr(Q) = Ai < 1, it holds that Ai = 1. As a result, we have 



Pr(Ci<i?) = Pr 



Pol 301 1 



ELiP^(^fc-l)-M|Elgfe||2 + l 
^ l + po||hu|p 



2« 



(35) 



Based on the result above, the problem of dllt is equivalent to the 
problem of il3\ . This completes the proof. 
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